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1 .  INTRODUCTION 


This  paper  discusses  the  use  of  pattern  recognition  to  automatically 
assess  targets  detected  by  the  installation  Security  Radar  (isr) 
system.*  The  ISR  is  a  uhf  site-security  radar  with  a  360-degree  elec¬ 
tronic  step-scan  capability.  A  coherent  receiver  develops  both  the  in- 
phase  and  quadrature  components  of  the  Doppler  return.  Detection  proc¬ 
essing  is  implemented  with  a  high-speed  microprogrammable  processor*  for 
1024  cells  chosen  from  the  available  8192  range -azimuth  cells.  Bach 
cell  is  11.5  degrees  wide  by  15  m  deep.  Target  classification  and  track 
processing  are  initiated  upon  target  detection. 

Radar  backscatter  from  a  range-azimuth  cell  contains  information 
about  the  scattering  objects  within  that  cell.  Target  classification  as 
addressed  in  this  pnper  uses  this  inherent  information  to  assign  priori¬ 
ties  to  possible  threat-target  tracks  and  to  minimize  the  formation  of 
tracks  on  nuisance  targets .  such  as  birds  and  animals,  and  false 
targets,  such  as  foliage  backscatter  and  system  noise.  Target  clas¬ 
sifiers  have  been  designed  and  evaluated  with  a  relatively  large  data 
base  obtained  jointly  by  Sandia  Laboratories  (SL)  and  the  Harry  Diamond 
Laboratories  ( HDL)  at  the  Department  of  Energy's  PANTEX  facility  in 
Texas.  The  data  were  obtained  with  a  Camp  Sentinel  Radar  (CSR)2  driving 
digital  recording  hardware  with  software  developed  at  SL. 


2,  PATTERN  RECOGNITION  APPROACH 

Conceptually,  pattern  recognition  involves  the  three  operations  of 
feature  definition  and  ranking,  feature  extraction,  and  class  discrimin¬ 
ation.  In  this  study  the  essential  target  features  were  determined  as 
derived  parameters  of  a  simple  radar-backscatter  return  model.  The 
Doppler  return  is  Fourier  transformed,  and  the  features  are  extracted  in 
the  spectral  domain.  Discriminant  functions  map  regions  of  the  N- 
dimensional  feature  subspace  to  regions  of  the  M-dimensional  decision 
space.  In  general,  various  classifier  structures  and  optimality 
criteria  are  available  for  selecting  the  class  discriminator.  in  this 
paper,  class  discrimination  is  based  on  the  value  of  the  conditional  a 
posteriori  probability  functions  evaluated  at  estimates  of  the  spectral 
features. 


•Concept  Evaluation  Phase  of  the  Installation  Security  Radar  System, 
HDL  internal  report,  24  December  1980  (available  from  author ). 

*R.  Lewis,  PMP-2  Processor  Handbook,  Stein  Associates,  Inc,,  Waltham, 
MA  (March  1978). 

2J.  Dent,  Camp  Sentinel  Radar  III  (U) ,  18th  Tri-Service  Radar 
Symposium  (1972).  (CONFIDENTIAL) 
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Although  numerous  distance  measures  have  been  suggested  for 
evaluating  feature  selection,  Kanal3  suggests  that  it  is  best  to  measure 
the  probability  of  error  directly.  For  the  isr,  the  a  priori  probabil¬ 
ity  of  each  target  class  is  not  known.  Therefore,  instead  of  maximising 
the  total  probability  of  error,  the  present  target  classifiers  were 
developed  by  maximising  the  probability  of  proper  classification,  while 
minimising  the  probability  of  improper  classification  for  each  target 
class  taken  separately. 

2.1  Target  Models 

The  use  of  physical  models  to  predict  the  spectral  character  of 
the  Doppler  return  offers  at  least  twrea  benefits  to  the  development  of 
target  classifiers.  First,  a  physical  model  allows  efficient  use  of  the 
available  data  base  since  only  the  essence  of  the  spectral  character¬ 
istics  of  the  return  is  used  to  develop  the  model  and  hence  design  the 
classifier.  Typically,  at  least  one-half  of  a  data  base  is  used  during 
classifier  training  (e.g.,  to  obtain  the  weights,  for  a  linear 
classifier).  Second,  strategies  to  determine  the  "best"  feature  subset 
are  not  necessary,  since  the  ranking  of  the  spectral  features  occurs 
naturally  from  the  development  of  the  target  model.  Third,  the  use  of 
target  modeling  yields  an  intuitive  relationship  between  the  spectral 
features  and  such  physical  parameters  as  target  velocity,  direction  of 
travel  relative  to  radial,  step  size,  and  step  rate.  Because  of  the 
structure  present  in  the  data,  this  type  of  target  modeling  (as  opposed 
to  less  intuitively  motivated  methods  such  as  principal-component 
modeling)  seems  to  be  well  suited  to  the  problem. 

2.1.1  Translating  Targets 

A  simple  three-point  target-backscatter  model  possesses  the 
essential  spectral  character  of  the  uhf  return  from  a  translating 
target.  The  modelled  return,  V,  is  the  vector  sum  of  the  contributions 
from  the  target  body  and  its  moving  appendages: 


3  4irr,  (t) 

V  -  E  Sj(t)  cos  - i -  .  (1) 

i-1  1  * 


The  a^(t)  (i  «  1,  2,  3)  represent  the  variation  in  cross  section  of  the 
scattering  point  as  a  function  of  time;  r^(t)  (i  *  1,  2,  3) 
represent  the  distance  between  the  ith  scattering  point  and  the  radar  as 
a  function  of  time;  and  \  is  the  rf  wavelength. 


3Ii.  Kanal  f  Patterns  in  Pattern  Recognition :  1963-1974 ,  IEEE  Trans. 

Inf.  Theory,  IT-20 ,  6  ( November  1974),  697-722 . 
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If  the  first  reflector  represents  the  return  from  the  torso, 
ri(t)  is  *  monotonic  function  of  time.  The  second  end  third  reflectors 
then  represent  the  return  from  the  moving  appendages  (e.g.,  human  arms 
and  legs,  bird  wings)  and  thus  r2<t)  and  r3(t)  possess  both  a  monotonic 
and  an  oscillatory  component.  The  return  from  a  wheeled  vehicle  can  be 
considered  a  degenerate  case  with  only  a  single  scattering  point.  For 
constant-velocity  targets,  V  can  be  expressed  as 

V  -  a  ;(t)  cos  (2flf,pt) 

+  a2<t)  cos  [2nf.pt  +  d>A( t )  ]  (2) 

+  a3(t)  cos  [2 irf.pt  4  4»2<t)  ]  , 


where 


2vr 

fT  -  — j*  cos  y 

VR  -  target  radial  velocity, 

y  «  angle  of  target  travel  relative  to  radial,  and 

<^(t)  and  $  (t)  are  angle  modulation  terms  caused  by  the 
moving  (resonant)  appendages. 


The  periods  of  a^ft)  (i  ■  1,  2,  3)  and  $u(t)  (j  ■  1,  2)  are  related  to 
the  appendage  resonant  frequency  (or  target  natural  frequency,  fo). 

Assuming  that  $^(t)  -  $  (t)  4  it  and  a2(t)  -  a3ft  +  V*o)» 
equation  (2)  can  be  written  in  the  simplified  form: 


a  ^(t)  cos  (2irf.pt  j  4  2a2(t)  cos  [2nf.pt  4  ^(t)  ]  .  (3) 


A  sample  of  the  measured  power  spectral  density  function  from 
a  radially  inbound  human  is  shown  in  figure  1 .  Figure  2  shows  an 
example  of  the  power  spectral  density  function  of  a  synthetically 
generated  target  signal  based  on  equation  (3).  The  modelled  human- 
target  return  (fig.  2)  can  be  seen  to  possess  much  of  the  same  spectral 
character  as  the  real  radar  return. 
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fHbUUENCY  (Hz) 

"T*0™  and  **"*  spectral  density 

100  ra?  J*9t  w*lkin?  radially  in-boundf 

steps  per  minute  with  normal  arm  and  leg  motion. 


function  of 

at  a  a  tap  rate  of 


FREQUENCY  (Hz) 

Figure  2.  Typical  Doppler  waveform  (no  range  cell  weighting)  and 
power  spectral  density  function  generated  by  generalized  target 
model  (eq  (3)). 


Human  return,— The  important  target  attributaa  that  appaar  to 
ba  obtainabla  from  tha  Dopplar  raturn  ara  (1)  targat  velocity, 
(2)  targat  natural  fraquancy,  (3)  distance  traveled  in  one  cycle  of  limb 
motion,  (4)  limb  modulation  characteristics,  (3)  indication  of  targat 
mass,  and  (6)  relative  sisa  of  tha  target*  Soma  of  these  attributes  can 
ba  estimated  direotly  from  tha  speotrum  of  tha  targat  raturn,  while 
others  must  ba  indirectly  estimated.  Because  of  shadowing,  roughness, 
multipath  effects,  the  target's  height  above  ground,  etc,  the  estimate 
of  relative  target  sise  may  be  in  error  by  more  than  an  order  of 
magnitude*  For  this  reason,  only  the  first  five  attributes  will  be 
consideredi 


radial  velocity,  VR>  i.e, 


target  velocity,  VT,  and  not 


where  y  is  typically  unknown.  A  possible  means  of  estimating  y  is  given 
shortly.  The  highest  velocity  for  a  sprinting  human  is  approximately  30 
ft/s  while  a  dog's  velocity  can  be  twice  that  high.  Thus,  although  a 
measured  velocity  of  greater  than  30  ft/s  does  not  inevitably  imply  a 
dog  target,  it  rules  out  the  possibility  that  the  return  is  due  to  a 
human . 


(1)  fT  is  related  to  the  target 


2V_ 


Of  considerably  mora  interest  is  the  actual 
just  its  radial  component >  i.e., 


VT  "  VR  *®e  Y 


(2)  fo  is  the  natural  frequency  of  the  target,  In  the  case 
of  a  human,  fo  is  equal  to  one  half  the  step  rate,  while  for  a  bird  it 
is  equal  to  the  wing-beat  frequency.  The  highest  expected  natural 
frequency  for  a  human  is  approximately  5  Hz,  while  a  duck  and  a  sparrow 
have  natural  frequencies  of  approximately  9  and  13  Hz,  respectively. 
Although  large  birds  have  wing-beat  rates  on  the  order  of  a  human's 
natural  frequency,  most  birds  have  significantly  higher  natural  frequen¬ 
cies.  In  the  case  of  nonnatural  targets  such  as  helicopters,  the 
propeller  rotation  rate  is  significantly  higher  than  most  bird  wing-beat 
rates. 


(3)  fT/fo  is  proportional  to  the  distance  of  radial  target 
travel  in  one  cycle  of  the  target's  natural  frequency.  This  distance 
measure  is  qiven  by  do,  where 


d 


0 


X  rT 

7  70- 


sec 


y  • 
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In  the  case  of  a  human,  do  represents  the  target  step  size.  An  estimate 
of  y  is  required  to  fully  use  the  quantity  fT/fo* 

(4)  The  nature  of  the  limb  modulation  is  expected  to  be  a 
species-dependent  characteristic#  In  addition,  it  has  been  empirically 
observed  that  by  comparing  the  relative  power  in  the  limb  modulation  and 
torso  return  for  natural  targets,  an  indication  of  the  angle  of  target 
travel  relative  to  radial  may  be  obtained.  Ihis  phenomenon  has  not  been 
fully  explained  but  may  be  due  to  masking  of  or  interference  among  the 
dominant  scatterers. 

For  a  human  on  a  radial  trajectory,  the  translation  lobe 
in  the  spectral  return  is  somewhat  larger  than  the  modulation  side- 
lobes.  As  the  angle  relative  to  radial  increases,  the  power  in  the 
translational  lobe  relative  to  the  power  in  the  modulation  harmonics 
diminishes  for  an  outbound  target  and  increases  for  an  inbound  target, 
as  shown  in  figures  3  and  4-,  respectively.  By  using  such  information, 
an  estimate  can  perhaps  be  made  of  the  direction  of  target  travel  rela¬ 
tive  to  radial.  An  estimate  of  y  could  be  used  by  both  the  tracking 
function  for  trajectory  prediction  and  in  the  calculation  of  target  step 
size. 


(5)  the  spectral  width  of  the  translational  component  and 
modulation  sidelobes  of  the  target  return  may  provide  additional  target 
discrimination.  As  suggested  by  L.  T.  James  (SL)  and  R.  Chase  ( HDL ) , 
this  spectral  wic  ,.h  contains  information  about  acceleration  of  the 
target  and  its  appendages.  In  certain  cases  the  spectral  width  may  be 
related  to  the  mass  of  the  target. 

For  an  upright  walking  human,  the  spectral  width  of  the 
target  translational  component,  AfT,  generally  increases  with  increasing 
target  velocity.  If  A fT  is  written  as  some  function,  g(),  of  the  target 
velocity,  VT,  then 


Assuming  that  the  inverse  of  g  exists  (i.e.,  g— 1 { ) ) ,  then  an  estimate  of 
Y  can  be  obtained  as 


Y  =  cos"1  — r  .  (5) 

Similarly,  since  do  =  (A/2)(fT/fo)  sec  y,  if  g-1()  exists,  do  can  be 
expressed  as 
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(6) 


Once  the  function  g  and  its  inverse  are  obtained,  the  estimates  of  y  and 
do  depend  only  on  the  easily  measured  spectral  attributes  fT,  fo,  and 
AfT.  The  two  separate  estimates  of  y  were  made  based  on  (1)  the  power 
in  the  translational  component  relative  to  the  modulation  sidelobes  and 
(2)  the  relative  width  of  the  translational  component.  These  single¬ 
cell  estimates  of  y  can  be  correlated  with  the  observed  target  trajec¬ 
tory  over  multiple  cells  to  produce  an  accurate  estimate  of  angle  of 
target  travel  relative  to  radial. 
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Figure  3.  UHF  return  from  an  outbound  human,  traveling  at  an  angle  of 
45°  relative  to  radial,  with  a  step  size  of  2  ft  and  step  rate  of  2 
steps  per  second. 
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Figure  4.  UHF  return  from  an  inbound  human,  traveling  at  an  angle  of 
45°  relative  to  radial,  with  a  step  size  of  2  ft  and  step  rate  of  2 
steps  per  second. 

Bird  return. — Birds  can  be  expected  to  present  a  problem  to  a 
ground-based  radar  designed  to  detect  humans  and  vehicles.  Tie  energy 
arriving  at  a  target  and  that  scattered  back  to  the  radar  consists  of 
both  a  direct  path  and  a  ground-reflected  wave.  A  vertical  lobe  struc¬ 
ture  in  the  radiation  pattern  results  from  the  interference  phenomenon, 
as  discussed  by  Skolnik,4  The  observable  effect  is  that  the  radar 
return  from  a  target  near  the  ground  falls  off  as  the  inverse  of  range 
to  the  eighth  power,  instead  of  the  inverse  of  range  to  the  fourth' 
power,  as  with  a  target  in  free  space. 


Skolnik,  Introduction  to  Radar  Systems,  McGraw-Hill  (1980),  Ch 
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For  path  length  differences  that  are  small  compared  to  a 
wavelength,  it  is  shown  in  appendix  A  that  the  received  power  from  a 
target  is  proportional  to 


f 


where 


ha  is  the  antenna  height, 

ht  is  the  target  height,  and 

R  is  the  (direct  path)  range  to  the  target. 

A  plot  of  the  approximate  variation  in  radar  return  power 
expected  from  a  point  target  as  a  function  of  target  height  is  shown  in 
figure  5.  For  example,  with  a  100-ft  high  radar  antenna,  the  return 
from  a  target  on  the  ground  would  be  approximately  12  dB  less  than  for 
the  same  target  at  a  height  of  100  ft. 

A  considerable  .’mount  of  work  has  been  reported  in  the  liter¬ 
ature  on  radar  return  from  bird  targets.  Birds  with  physical  extent 
somewhat  less  than  one-half  of  the  rf  wavelength  appear  as  Rayleigh 
(-catterers  with  cross  sections  exhibiting  a  X”*4  wavelength  dependence. 
At  wavelengths  less  than  the  target  extent,  the  cross  section  is 
expected  to  be  proportional  to  the  wavelength.5  At  approximately  3  GHz, 
the  cross  section  of  all  birds,  regardless  of  size,  exhibits  minimum 
spread  (resonance).5  Figure  6  compares  the  mean  cross  section  of  a 
human  to  that  of  several  birds  in  the  range  of  0.4  to  10  GHz.  Table  1 
lists  available  average  cross-section  measurements  at  uhf  (approximately 
400  MHz)  for  several  species  of  birds.  Most  sources  suggest  that  the 
amplitude  of  the  bird  return  follows  a  lognormal  probability  density 
function.6  Cross-section  measurements  of  birds  in  flight  performed 
simultaneously  at  uhf,  S,  and  X  band  have  shown  the  standard  deviation 
of  the  return  to  be  of  the  same  order  of  magnitude  as  the  mean  of  the 
return . 6 


The  most  important  discriminant  reported  in  bird  return  is  an 
amplitude  modulation  whose  fundamental,  fo,  is  the  wing-beat 


%.  Pollon,  Distribution  of  Radar  Angels,  IEEE  Trans .  Aerosp. 
Electron.  Syst.,  AFS-8 ,  6  (November  1972),  721-727. 

Atlas  and  K .  Hardy,  Radar  Analysis  of  the  Clear  Atmosphere: 
Angels,  Progress  in  Radio  Science  1963-66,  International  Scientific 
Radio  Union,  15th  General  Assembly,  Munich  (5-15  September  1966). 


16 


4i  -Kfcr. 4i>* > i.ililMkuii dM. ****.1*. n,  ..,<11  «iUiJhiM<4**ll«U «lci 


MEAN  RADAR  CROSS  SECTION  (dBsm) 


FREQUENCY  (GHz) 

Figure  6.  Examples  of  radar  cross  section  of  birds  (excerpted  from  ref 


T 

.4 

s 


Wing-beat  frequencies  of  between  11  end  15  beets  per  second 
were  measured  for  e  number  of  flights  of  e  smell  pheesent  et  the  PANTEX 
facility.  The  results  egree  reesonebly  well  with  the  published  results 
for  comparable-sise  birds  shown  in  teble  3. 


TABLE  1.  UHF  BIRD  RADAR  CROSS  SECTION  (ref  7) 


Target 

Reder 

cross  section 
(dBsm) 

Source 

Pigeon 

-30 

Technology  Service  Corp.* 

-32 

Condon  ^ 

Grackle 

-42 

Technology  Service  Corp. 

-44 

Condon 

Sparrow 

-55 

Technology  Service  Corp. 

-57 

Technology  Service  Corp. 

aA  Clutter  Model  for  Artillery  and  Mortar  Locatiirj 
Radar,  Technology  Service  Corp .,  TSC-PD-041-1 ,  prepared 
for  Nitre  Corp .  and  USASCOM  (7  July  1969). 

bE.  J.  Condon,  Scientific  Study  of  Unidentified 
Flying  Objects ,  New  for*  Tims,  Bantam  Books  (1968), 
655-716. 


TABLE  2.  SUGGESTED  MODELS  RELATING 


WING  LENGTH  AND  CHARACTERISTIC 

WING-BEAT  RATE  (ref  9) 

Model* 

source 

f  0 

1  1.15  «  3540 

Greenewalt  b 

fo 

t  9.827  -  572 

Houghton  ° 

af0  is  the  wingbeat 

rate  (Hz)  and  l  is  the 

length  of  the  wing  (cm). 

“C.  H.  Greenewalt,  Hummingbirds,  Doubleday 
i  and  Co.,  Inc .,  Garden  City,  NY  (1960). 

*  CE.  W.  Houghton,  The  .Effect  of  Changes  in 

^  Target  Trajectory  on  the  Wlngbeat  Modulation 

J  Pattern,  Royal  Radar  Establishment,  RRE  Memo- 

•  randum  no.  2456,  Malvern,  Worcs .,  England  (1969). 


Sine#  individual  birds  or  TABLE  3.  TYPICAL  WING-BEAT 
flocks  of  birds  could  pr ^sent  serious  BATES  (raf  7) 

false-alarm  problems  to  a  security  radar  — . —  ■■  —  —  -  — — 

system  (especially  at  times  of  high 

activity  such  as  the  early  mortiing,  late  Species  Rate  (beats/s) 
evening,  and  migration  seasons),  eonsid-  ■  -  — — 


erably  more  backscatter  return  from  bird 

Sparrow 

13 

targets  needs  to  be  obtained  and 
analysed  in  the  future.  Depending  on 

Duck 

9 

the  size  of  the  bird  and  whether  it  is 
taking  off,  landing,  flying  level,  or 

Crow 

3  to  4 

soaring,  various  wing-beat  rates  can  be 
expected.  Birds  generally  change  direc¬ 

Stork 

2 

tion  more  often  and  have  a  higher  trans¬ 
lational  velocity  (fig.  7)  and  a  higher 

Pelican 

1.1 

natural  frequency  (table  3)  than  a  human 
target;  thus,  by  monitoring  the  change 

Hummi ngbl rd 

10  to  80 

in  fT#  fo,  and  fT/fo  along  a  target  ■■■  - . ■■■■ .  . 

trajectory  (through  one  or  more  range-azimuth  cells),  reliable 
algorithms  can  perhaps  be  developed  to  classify  those  tracks  caused  by 
birds. 


Animal  return.— Target  models  similar  tr  the  human-target 
model  need  to  be  developed  for  each  animal  specie-  that  is  expected  to 
cause  nuisance  alarms  at  a  given  radar  site.  Since  the  uhf  return  from 
a  biped  can  apparently  be  adequately  modelled  as  an  angle-modulated 
return  with  varying  intensity  depending  on  the  step  rate,  quadrupeds 
should  produce  an  intensity  and  angle  modulation  with  a  fundamental  that 
corresponds  to  one-half  or  one-fourth  the  animal's  step  rate. 

Vehicle  and  aircraft  return. — Wheeled  vehicles  such  as  cars 
and  small  trucks  appear  to  give  an  essentially  monotone  return,  whereas 
tracked  vehiclea  may  produce  return  that  contains  certain  modulation 
harmonics.  Helicopters  and  propeller-driven  aircraft  are  expected  ■  »> 
produce  return  with  modulation  harmonics  related  to  the  propeller  rota¬ 
tion  rate. 


Insect  return. — Swarms  of  flying  insects  have  long  been 
suspected  as  one  of  the  causes  of  radar  angels  (echoes  observed  from 
apparently  clear,  targetless  regions).  Attempts  to  measure  individual 
insect  radar  cross  sections  and  spectral  characteristics  of  their  return 
have  been  generally  carried  out  at  relatively  short  wavelengths. 
However,  at  least  one  set  of  cross-section  measurements  was  performed  at 
both  uhf  and  L  band.10  Table  4  gives  empirically  determined  target 

l(V,  Glover,  R,  Hardy,  T .  Konrad,  W .  W.  Sullivan,  and  A.  Michaels, 
Radar  Observations  of  Insects  in  Free  Flight,  Science,  1 54 ,  3752  (25 
November  1966),  967-972, 
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1 


u 

a 


croaa  aactiona  of  thraa  common  Intacta  at  L  band*  while  at  uhf  (inter¬ 
estingly  enough  at  435  MHa— approximately  the  same  frequency  as  the  CuR) 
the  target  eroaa  aactiona  ware  below  the  minimum  detectable  cross 
section*  For  wavelengths  greater  than  10  cm*  an  upper  bound  on  the 
insect  cross  section  ia  established  by  a  X"3  relationship. 10  The  insect 
croas  sections  appearing  in  table  4  would  be  diminished  by  at  loaat  an 
order  of  magnitude  at  433  MHs. 
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Figure  7.  Range  of  bird  velocities  (abstracted  from  ref  7). 


TABLE  4.  L-BAND  RADAR  CROSS  SECTION  rr 
THREE  INSECTS  MEASURED  T”  LREE 


FLIGHT 

(ref  10) 

Insect  target 

L-band  radar 
cross  section  (dBsm) 

wingless  Hawkmoth 

-60 

Honeybee 

-63 

Dragonfly 

-67 

Since  no  studies  have  been  made  to  determine  the  statistics 
of  insect  return,7  there  is  little  means  of  determining  whether  high 
enough  insect  densities  will  be  found  at  a  given  geographical  location 
to  cause  a  detection  and  subsequent  track  formation. 

Because  there  is  motion  of  the  insects  within  a  swarm,  the 
spectral  return  may  not  provide  a  good  measure  of  the  swarm's 
velocity.  However,  if  the  fundamental  of  the  wing-beat  rate  is 
observed,  it  could  be  used  to  discriminate  bird  from  insect  return  since 
insect  wing-beat  rates  arc  aenerally  higher  than  bird  wing-beat  rates. 
For  example,  locust  wing-beat  frequencies  range  from  14  to  20  Hz, 
butterflies  from  20  to  30  Hz,  and  crows  from  only  3  to  4  Hz. 

2.1.2  Nontranslating  "Target"  Return 

The  power  spectral  density  function  of  the  complex  (quad¬ 
rature)  return  due  to  windblown  objects  such  as  foliage,  power  lines, 
etc,  that  make  no  net  progress  toward  the  radar,  consists  of  approxi¬ 
mately  equal  positive-frequency  and  negative-frequency  power%  in  most 
cases  the  power  spectral  density  function  of  the  return  is  nearly  sym¬ 
metrical  about  0  Hz. 

An  empirically  determined  model  that  predicts  the  long-term 
average  power  spectral  density  function  for  uhf  and  L-band  windblown- 
foliage  backscatter  return  was  previously  developed  at  HDL,11  A  simple 
exponential  was  found  to  describe  the  data  for  wind  conditions  between  2 
and  20  mi /hr.  The  expression  for  the  power  spectral  density,  P(f),  is 

7 ? Reedy  and  T.  Cutler,  the  HOWLS  Radar  Sky-Clutter  Environment , 
Engineering  Experiment  Station,  Georgia  Institute  of  Technology,  ESD-TR - 
75-319  (15  September  1975). 

AiR.  Antony,  C.  Roberts,  and  s.  Peperone,  An  Investigation  of  Radar 
Backscatter  from  Foliage ,  18th  Tri-Service  Radar  Symposium  (1972). 
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2.2  Target  Features 


The  important  target  characteristics  apparent  from  the  trans* 
lation  target  model  are  the  measured  radial  frequency,  fT,  and  the 
claas-characteriatic  amplitude  modulation,  a^(t),  and  angle  modulation, 
^(t).  The  following  eix  features  have  been  derived  from  equation  (S) 
and  were  discussed  in  the  last  section  as  applicable  to  the  various 
target  models > 


(a) 

fT 

(b) 

fo 

(c) 

do 

(d) 

nS 

(e) 

Af  »p 

(f) 

Afo 

target  radial  translational  frequency 

target  natural  frequency  (e.g,,  half  the  step  rate, 
RO,  of  a  human;  wing-beat  rate  of  a  bird) 

(AVdfo)  sec  y  (e.g.,  step  size  of  human,  distance 
of  travel  in  one  wing  beat  of  a  bird) 

number  of  harmonics  of  fo  observed  in  the  power 
spectral  density  function  of  the  return  within 
specified  power-level  and  signal-to-noise  ratio 
constraints 

spectral  width  of  the  target  translational 
frequency 

spectral  width  of  the  fundamental  of  the  target 
resonant  frequency 


The  feature  ns  accounts  for  the  gross  nature  of  the  class-characteristic 
amplitude  and  angle  modulation  observed  in  the  return.  The  feature  AfT 
is  related  to  the  rate  of  change  of  the  target  translational  velocity 
for  straight-line  motion.  Target  acceleration  is  generally  greater  for 
a  low-maos  target  such  as  a  bird  than  for  a  relatively  high-mass  target 
such  as  a  human.  The  spectral  spreading  of  the  target  resonant 
frequency  tends  to  be  a  class -dependent  characteristic  that  is  useful  in 
high  signal-to-noise  ratio  return. 


The  values  of  features  are  estimated  by  algorithms  that 
associate  the  features  with  the  20  largest  distinct  peaks  of  the  power 
spectral  density  function  of  the  return.  Proper  classifier  performance 
depends  on  the  correct  determination  of  f rp  and  fo*  Since  a  target  can 
be  incoming  or  outgoing,  the  spectral  character  of  the  nontranslating 
part  of  the  target  return  is  biased  by  the  radial  translational 
frequency  so  that  most,  but  not  all,  of  the  return  consists  of  either 
positive  or  negative  frequency  components.  Since  the  limb  motion  occurs 
forward  and  backward  of  the  torso,  the  return  should  have  nearly  equal 
power  in  both  the  upper  and  lower  sidebands  of  the  translational  lobe. 
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This  criterion  has  been  used  effectively  to  select  the  proper  value  of 
fT  even  in  cases  where  some  spectral  sidelobes  are  an  order  of  magnitude 
larger  than  the  translational  lobe. 

The  target  velocity  depends  on  the  product  of  the  target  step 
sice,  do»  and  the  target  step  rata,  Rq.  Assuming  that  the  target  step 
9isa  is  directly  proportional  to  the  target  step  rate  (the  faster  one 
walks,  the  larger  the  step  sice),  then  a  "consistency  criterion"  can  be 
developed  to  associate  the  most  likely  spectral  peak  with  the  spectral 
feature  fo. 

2.3  Target  Discrimination 

Target  discrimination  as  addressed  in  this  paper  is  based  on 
the  classical  likelihood  ratio  for  composite  hypotheses  (parameters  that 
are  random  variables)! 


J  p(y|eaHA)  pMhJ  dp 

ij  /  pCylSfHj)  pfolfij)  d©  1 

where  0  is  a  vector  random  variable  representing  the  set  of  target 
features  as  discussed  in  section  2.2,  and  designates  target  class 
i.  When  an  estimate  of  the  feature  vector,  0,  has  been  made,  the  like¬ 
lihood  ratio  reduces  to 


pfelHj) 

ij  "  pfyje.Hj)  ptolHj) 


(9) 


Applying  Bayes*  rule,  equation  (9)  can  be  rewritten  in  the  form 


A 


ij 


p(y|Hj) 

p(©|y#Hj)  p(y|Hj  ) 


p(q1y>hi) 

p(o |y»  «j ) 


AijCyl®) 


(10) 

(11) 


where  A^a(y|0)  is  the  simple  hypothesis  (fixed  parameter)  likelihood 
ratio  ana  p(o|y,HjJ  is  the  a  posteriori  probability  density  function 
(pdf)  of  the  parameter  vector  conditioned  on  the  target  class,  evaluated 
at  the  estimates  of  the  target  features. 
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The  objective  ie  to  develop  the  a  posteriori  pdf  to  be 
evaluated  at  0  ■  Si  The  elements  of  0  are  defined  as  the  set  of 
spectral  features  described  in  the  last  seotiont 

©  ■  (f^i»f 0»ng»df^i» Af o)  • 

For  tractability,  the  joint  a  posteriori  pdf  will  be  approximated  by 

p(o|y»Hi)  -  p(0i|y»Ht)  p(©2|y»H1)  p(0a|y,Hi)  ,  (u) 


where 

A  /  > 

0i  ■  ( f  ij>,  f  0 » ng  J , 

A  t  ^ 

©2  “  (  f  IJI»  AfijJ  , 

A  .  , 

03  ■  (fo,AfoJ ,  and 

pf  01  I  y  #  HA )  -  p(fT|  Ha)  pf  io  |  Ha  )  P^oIhJ  pfbglHj  pfP.jHj  . 

The  univariate  a  priori  pdf 'a  of  the  features  were  developed  based  on 
physical  constraints,  published  data,  and  empirical  measurements.  The  a 
priori  probability  p(do|H^)  explicitly  expresses  the  dependency  between 
target  velocity  (i.e.,  fT)  and  target  resonant  frequency  (i.e.,  f0). 
The  pdf  P (pan I Hi  )  relates  to  the  confidence  in  the  spectral  feature 
estimates  based  on  the  measured  signal-to-noise  ratio  of  the  return. 

In  a  similar  manner, 

p(©2|y.Hi) 

and 

p(03|y,Hi) 


As  an  example,  the  probability  functions  of  the  seven  univar¬ 
iate  feature  density  functions  for  a  human  target  are  presented  in 
figure  9,  Similar  conditional  pdf's  apply  to  all  other  translating 
target  return. 

The  class-discriminator  decision  rule  is,  then,  assign  the 
return  as  due  to  that  class  with  the  largest  value  of  p(©|y,H^).  In 
general,  this  represents  a  nonunique  assignment  to  a  target  class,  since 
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P (o| y^Ht  j  could  be  the  same  for  more  than  one  class.  An  advantage  of 
this  decision  rule  is  that  it  quantifies  the  likelihood  of  proper  clas¬ 
sification  and  indentifies  other  target  classes  possessing  nonzero 
probabilities. 


SNR 


Figure  9,  Univariate  probability  density  functions  for  a  human 
target. 


3.  CLASSIFICATION  RESULTS 


3.1  Single-Cell  Assessment 

The  extent  of  the  relatively  high  signal-to-interference  ratio 
target  data  base  used  for  classifier  design  and  testing  is  summarized  in 
table  5.  The  approximately  330  files  of  human-target  return  include 
data  for  both  single  and  multiple  humans.  The  single-human  return 
includes  walking,  jogging,  and  running  targets  with  various  step  sizes, 
step  rates,  and  aspect  angles.  The  multiple -human  return  includes 
targets  walking  in  step  and  out  of  step,  with  some  carrying  Ml  6 
rifles.  A  medium-sized  dog  (40  lb)  was  used  to  obtain  the  only  animal 
data  presently  available.  A  pigeon-sized  pheasant  was  released  from 
within  a  range-azimuth  cell  to  obtain  the  bird  backscatter  data. 

Examples  of  the  power  spectral  density  functions  of  the  target 
return  from  a  radially  inbound  walking  and  running  human  are  shown  in 
figures  10  and  11,  respectively.  An  example  of  an  inbound  human  walking 
at  an  angle  of  45  degrees  relative  to  radial  is  shown  in  figure  12.  The 
power  spectral  density  functions  for  three  inbound  humans,  an  inbound 
dog,  and  an  inbound  vehicle  are  shown  in  figures  13  through  15, 
respectively.  Finally,  an  example  of  the  power  spectral  density 
function  of  the  backscatter  from  wind-blown  power  lines  is  shown  in 
figure  16. 


TABLE  5.  SUMMARY  OF  ORIGINAL  DATA  BASE 


Target 

Number  of 
experiments 

Single  man 

(walking,  jogging,  running) 

306 

Multiple  man 
(walking,  running) 

23 

Dog 

26 

Small  pheasant 

9 

Clutter 

50 

Vehicles 

19 

28 


TIME  (s| 


CLASSIFIER  OUTPUT. 
PROS  CLUTTER  =  0.00 
PROB  VEHICLE  =  0.00 
PROB  HUMAN  =  1.00 
PROB  ANIMAL  =  0.00 
PROB  BIRD  =  0.00 

SPECTRAL  FEATURES: 

f,  -4.19 

*  0  '  194 
d0  *  2-6 


FREQUENCY  (Hz) 


Figure  10.  Radar  backscatter  return  from  an  inbound  walking  human  and 
corresponding  power  spectral  density  function;  return  was  classified  as 
due  to  a  human  with  probability  1.0. 


The  performance  of  a  classifier  can  be  expressed  in  terms  of 
the  matrix  where  pfH^|Hj)  is  the  probability  of  choosing 

hypothesis  H^,  given  that  Hj  is  true.  If  this  matrix  is  diagonal,  the 
classification  would  be  ideal.  The  results  of  the  single-cell 
assessment  using  the  current  target  classifiers  for  the  present  data 
base  are  shown  in  table  6.  For  the  compilation  of  these  data,  the 
return  was  classified  as  due  to  that  target  class,  H^,  whose  probability 
assessment,  p(©|y,HjJ,  was  greater  than  or  equal  to  0.5. 
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Figure  11.  Radar  backscatter  from  an  inbound  running  human  and 
corresponding  power  spectral  density  function;  return  was  classified  as 
due  to  a  human  with  probability  1.0. 
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Figure  12.  Radar  backscatter  from  an  inbound  human,  walking  at  a  45° 
angle  relative  to  radial,  aud  corresponding  power  spectral  density* 
return  was  classified  as  due  to  a  human  target  with  probability  1.0. 
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Figure  13*  Radar  backscatter  return  from  three  inbound  humans  and 
corresponding  power  spectral  density  function;  return  was  correctly 
classified  as  due  to  a  human  with  probability  1*0. 
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Figure  14.  Radar  backscatter  from  a  40-lb  dog  and  corresponding  power 
spectral  density  function;  return  was  classified  as  due  to  an  animal 
with  probability  1.0  and  as  a  human  with  probability  0.01. 
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Figure  15.  Radar  backscatter  from  an  Inbound  car  and  corresponding 
power  spectral  density  function;  return  was  classified  as  a  vehicle 
wi th  probabi li ty  1.0. 
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Figure  16.  Radar  backscatter  from  windblown  power  lines  and 
corresponding  power  spectral  density  function;  return  was  classified  as 
due  to  nontranslating  target  with  probability  1.0. 


TABLE  6.  SINGLE-CELL  CLASSIFIER  PERFORMANCE  SUMMARY 


Classification 


Actual 

target 

Clutter 

Vehicle 

Human 

Animal 

Bird 

Clutter 

1 .0 

0 

0 

0 

0 

Vehicle 

0 

0.95 

0 

0.05 

0 

Single  man 

0 

0 

0.92 

0.08 

0 

Multiple  men 

0 

0 

0.78 

0 

0 

Dog 

0 

0 

0.04 

0.88 

0.11 

Pheasant 

0 

0 

0 

0.11 

0.89 

Bor  the  available  data,  no  nontranslating  "target"  return  was 
misclassified,  nor  was  any  translating  target  return  classified  as  due 
to  nontranslating  return.  The  probability  of  a  human-target  return 
being  properly  classified  was  approximately  92  percent*  most  of  the 
remaining  8  percent  of  the  human-target  return  that  was  not  properly 
classified  was  from  jogging,  sprinting,  and  crouching  targets.  Approx¬ 
imately  4  percent  of  the  dog  return  was  improperly  classified  as  due  to 
a  human.  One  pheasant  return  was  improperly  classified  as  due  to  an 
animal.  Nearly  89  percent  of  the  dog  return  was  classified  as  due  to  an 
animal.  About  5  percent  of  the  vehicle  return  and  8  percent  of  the 
human  return  was  improperly  classified  as  due  to  an  animal. 


Multiple-Cell  Assessment 


Although  single-cell  target  classification  may  not  be  entirely 
adequate,  as  shown  in  table  6,  track  classification  based  on  the  full 
track  history  offers  a  means  of  reducing  the  probability  of  incorrect 
target  classification. 


One  possible  approach  to  designing  a  track  classifier  is  to 
assume  that  the  single-cell  assessments  are  independent  from  cell  to 
cell  within  the  track.  The  statistics  collected  for  the  single-cell 
assessment  involve  the  same  targets  making  repeated  passes  through  a 
single  range  gate.  The  extension  from  multiple  passes  through  a  single 
cell  to  a  single  pass  through  multiple  cells  is  thus  not  entirely  unrea¬ 
sonable  . 
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Assuming  independence  from  cell  to  cell,  table  7  tabulates  the 
probability  of  at  least  L  classifications  as  a  human  target  for  a  track 
of  length  K  cells.  For  instance,  approximately  4  percent  of  the  dog 
return  was  classified  as  due  to  a  human  target;  however,  the  probability 
that  a  target  will  be  classified  as  a  human  at  least  2  out  of  3  times  is 
98  percent  if  the  target  is,  in  fact,  a  human  and  essentially  sero  if 
the  target  is  a  dog.  Similar  multiple-cell  assessments  for  classifi¬ 
cation  as  an  animal  or  a  bird  are  given  in  tables  8  and  9.  Thus,  based 
on  the  results  of  the  current  classifiers  operating  on  the  present  data 
base,  high-reliability  target  assessment  for  tracks  as  short  as  3  or  4 
contiguous  cells  may  be  possible. 


TABLE  7.  PROBABILITY  OF  AT  LEAST  L  CLASSI¬ 
FICATIONS  AS  HUMAN,  GIVEN  K  TOTAL 
CELLS  ASSESSED 


Actual 

target 

L  -  1 

K  •  1 

L  -  1 

K  «  2 

L  -  2 

K  -  3 

Human 

0.92 

0.99 

0.98 

Dog 

0.04 

0.08 

0.00 

TABLE  8, 

PROBABILITY  OF  AT  LEAST  L  CLASSI¬ 
FICATIONS  AS  ANIMAL,  GIVEN  K  TOTAL 
CELLS  ASSESSED 

Actual 

L  «  1 

L  -  1 

L  -  2 

target 

K  ■  1 

K  «  2 

K  -  3 

Dog 

0.88 

0.92 

0.96 

Pheasant 

0.11 

0.01 

0.00 

Human 

0.08 

0.00 

0.00 

Vehicle 


0.05 


0.00 


0.00 


TABLE  9.  PROBABILITY  OF  AT  LEAST  L  CLASSIFICATIONS 
AS  BIRD,  GIVEN  X  TOTAL  CELLS  ASSESSED 


Actual 

target 

L  -  1 

X  -  1 

L  -  1 

X  -  2 

L  -  2 

K  •  3 

L  -  3 

X  -  4 

Bird 

0.89 

0.99 

0.96 

0.94 

Dog 

0.11 

0.21 

0.03 

0.00 

A  single  target  may  be  detected  in  multiple  range -azimuth  cells 
because  of  range  sidelobes,  partially  overlapping  azimuthal  antenna 
patterns,  or  multipath  effects*  Although  operations  such  as  simple 
centroiding  can  be  used  to  reduce  these  multiple  detections  to  a  single 
range-azimuth  cell  to  be  added  to  the  track  vector,  signature  analysis 
of  the  return  can  be  used  for  this  purpose.  Since  the  suggested  sig¬ 
nature  processing  isolates  the  target  radial  velocity  component  in  the 
spectrum  of  the  return,  the  probable  target  location  can  be  assigned  to 
that  cell  with  the  largest  radial  translational  component.  In  addition, 
the  separation  of  multiple  targets  with  different  radial  velocities  and 
return  with  clearly  different  spectral  character  is  possible  before 
simple  centroiding  reduces  multiple  alarming  cells  to  a  single  alarm. 

A  historical  assessment  of  the  spectral  features  of  the  full 
track  history  may  also  provide  additional  discrimination.  A  case  in 
point  would  be  for  a  bird'  target.  While  a  bird  may  soar  through  any 
given  range-azimuth  cell  on  a  path  nearly  tangential  to  the  radar,  it 
may  subsequently  change  to  a  more  radial  path  and  begin  beating  its 
wings.  Thus,  although  the  initial  assessment  of  the  spectrum  of  the 
return  would  yield  a  poor  estimate  of  the  target's  actual  velocity  and 
no  apparent  modulation  at  the  species'  characteristic  wing-beat  rate, 
subsequently  the  return  might  be  completely  characteristic  of  bird 
return . 


Thus,  the  overall  track  assessment  could  be  based  on  the  inde¬ 
pendent  classification  assessments  from  cell  to  cell  within  the  track 
and  on  the  statistics  of  the  measured  spectral  features. 


Tha  track  assessment,  Ptrack'  for  K  total  calls  assessed  night 
ba  of  tha  form 


ptraek^ 


Bo(m) 

K 


K 

l 

k-1 


Pfctm) 


N 


+  E 
n»1 


^(m)  py(n,m) 


03) 


+  $,(m)  p0(n,m)  , 


whare 

m  is  tha  target  class  (M  total  classes); 

pj, (m )  is  the  probability  assessment  of  tha  kth  cell  within  the 
track  (k  -  1 ,  . . . ,  K) j 

p^O  is  the  probability  assessment  for  the  mean  of  the  nth 
spectral  feature  for  each  target  class  (n  ■  1 ,  .  .  . ,  N)i 

pg()  is  the  probability  assessment  for  tha  standard  deviation 
of  the  nfc"  spectral  feature  for  each  target  class;  and 

Bo*  Bj^*  are  constants  that  depend  on  the  target  class  m. 

The  track  is  classified  as  that  target  class  m  with  the  largest  value  of 
ptrack^  * 


For  a  bird  target  under  track,  although  several  p^fbird)  might 
be  nearly  aero  if  the  bird  were  soaring  through  these  cells,  a  high 
standard  deviation  for  the  spectral  feature  fo  for  the  entire  track,  for 
instance,  would  be  assigned  a  high  value  of  pa(bird).  the  track  would 
be  classified  as  due  to  a  bird  target  if 

« 

ptrack^r<s)  £  ptrack^m^ 

for  all  m  with  associated  confidence  Ptrack^*^  *  it  would  be  clas¬ 
sified  as  an  unknown  target  if  all  Ptracfcrm)  <  Pmin'  where  pmin  is  a 
constant.  Additional  feature  statistics  such  as  the  minimum  and  the 
maximum  of  a  given  feature  over  the  entire  track  history  may  provide 
further  target  discrimination. 

An  example  of  the  power  spectral  density  function  of  the  return 
from  a  wild  duck  flying  through  three  consecutive  range-azimuth  cells  is 
shown  in  figure  17.  The  bird  beats  itB  wings  while  in  range  cell  No.  1 
flying  at  a  radial  speed  of  42  ft/s,  sets  its  wings  in  range  cell  No.  2, 
and  changes  direction  and  soars  out  of  range  cell  No.  3.  This  pattern 
is  apparent  from  both  the  time  series  and  power  spectral  density 
function  of  the  returns.  On  a  single-cell  basis,  the  bird  was  correctly 
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classifiad  in  ranga  call  No.  1,  misolassifiad  as  a  vahicla  in  aall  No,  2 
(baeausa  of  tha  lack  of  a  wing-baat  modulation),  and  claaaifiad  as  an 
unknown  targat  in  call  No,  3*  Though  on  a  aingla-call  basis,  tha  raturn 
is  miaclassifiad  in  ranga  calls  No,  2  and  3,  tha  targat  track  might 
still  ba  proparly  classifiad  sinca  such  a  pattarn  is  charactaristic  of 
birds  in  flight. 
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Figure  17  (b).  Return  was  classified  as  due  to  a  vehicle  because  of 
lack  of  wing-beat  harmonics  in  cell  2. 


41 


4.  INTEGRATED  DETECTION/CLASSIFICATION  SYSTEM  SIMULATION 

A  digital  simulation  of  the  detection  processor  has  been  developed 
that  allows  the  processing  of  real  radar  return  on  a  gene rel -purpose 
computer.  This  permits  the  development  of  data  handling  and  processing 
software  which  will  not  interfere  with  the  radar  hardware.  A  block 
diagram  of  the  system  interaction  is  shown  in  figure  18.  The  outputs  at 
all  stages  of  processing  can  be  monitored  on  a  graphics  display. 


Figure  18.  Block  diagram  of  proposed  radar  system. 


when  a  cell  alarm  occurs,  the  raw  Doppler  return  from  that  cell  is 
collected  for  the  estimated  time  the  target  will  remain  within  the 
range-azimuth  cell.  An  estimate  of  this  interval,  t,  is  selected  as 

»  2aAR 

T  =  -  ,  (14) 

fRl 

where  AR  is  the  range  cell  depth,  fR  is  the  estimated  radial  trans¬ 
lational  frequency,  the  parameter  a  is  a  dimensionless  quantity  that 
depends  on  the  signal-to-noise  ratio  and  various  processor  parameters, 
and  A  is  the  rf  wavelength.  An  estimate  of  fR  is  obtained  from  the  fast 
Fourier  transform  of  the  first  128  samples  of  the  collected  data. 

A 

Once  t  seconds  of  raw  Doppler  from  an  alarming  range-azimuth  cell 
have  been  collected,  the  data  are  weighted  with  a  Kaiser-Bessel  window 
and  a  512-point  discrete  Fourier  transform  (DPT)  is  computed.  The 
magnitudes  of  the  complex  spectral  coefficients  are  calculated  and  the 
20  largest  distinct  peaks  of  the  resulting  estimated  power  spectral 
density  function  are  determined.  The  spectral  lobes  are  then  associated 
with  the  various  spectral  features  and  the  target  classifiers  produce 
probability  assessments  for  each  of  the  five  target  classes. 


FILTER  NUMBER 


An  example  of  the  processor  output  for  a  dog  passing  through  a 
single  range-azimuth  cell  is  shown  in  figure  19.  In  this  case,  5  of  the 
12  balanced  processor  channels  were  monitored.  The  fixed  threshold 
level  was  purposely  set  high  so  that  no  threshold  crossing  occurred. 


TIME  {•) 

Figure  19.  Five  of  first  seven  simulated  high-velocity  filter  outputs 
for  a  dog  target,  passing  outbound  through  a  single  range-azimuth  cell; 
fixed  threshold  was  set  so  that  a  threshold  crossing  would  not  occur. 


In  figure  20,  the  same  data  are  processed  using  an  adaptive  constant 
false-alarm  rate  threshold  level.  In  this  case,  detection  occurs  and 
the  incoming  data  begin  filling  a  temporary  buffer.  A  128-goint  DFT  is 
calculated  and  a  refined  estimate  of  t  is  computed.  After  t  seconds  of 
the  Doppler  return  has  been  collected,  the  target  classification 
algorithms  are  initiated.  Figure  21  shows  the  Doppler  waveform  of 
length  t  seconds,  the  corresponding  estimate  of  its  power  spectral 
density  function,  and  the  classification  results.  In  this  case,  the 
detection  was  assessed  as  due  to  an  animal  target  that  is  outbound  with 
a  radial  translational  frequency  of  22.34  Hz  and  a  target  natural 
frequency  of  3.27  Hz. 
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Figure  20.  Simulated  high-velocity  filter  outputs  for  a  dog  target, 
passing  outbound  through  a  single  range-azimuth  cell;  adaptive 
threshold  was  used. 


Similar  results  for  both  human  and  bird  targets  are  shown  in  figures 
22  to  26.  In  all  cases,  the  target  classifiers  correctly  assessed  the 
cause  of  the  processor  alarms. 

Finally,  the  processor  output  for  the  return  from  a  range-azimuth 
cell  containing  no  target  is  shown  in  figure  27.  In  this  case,  the 
binary  integrator  output  does  not  cross  thr  threshold  and  consequently 
no  target  detection  occurs. 
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Figure  21.  Radar  backscatter  return  from  a  running  dog,  collected 
following  target  detection  indie. '  1  in  fig.  20>  estimated  duration  of 
target  in  range-azimuth  cell,  t,  v;as  calculated  to  be  approximately 
1.8  s . 


TIME  |o) 

Figure  22.  Simulated  low-velocity  filter  outputs  for  a  human  passing 
outbound  through  a  single  range-azimuth  cell)  fixed  threshold  was  set 
so  that  a  threshold  crossing  would  not  occur. 
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Figure  25.  Simulated  high-velocity  filter  outputs  for  a  pheasant  flying 
inbound  through  a  single  range-asimuth  celli  adaptive  threshold  was 
used. 
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Figure  26.  Radar  backscatter  return  from  a  pheasant,  collected 
following  target  detection  indicated  in  fig.  25i  estimated  duration  of 


target  in  range-asimuth  cell  was  calculated  to  be  approximately  0.9  s. 
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Figure  27.  simulated  high-velocity  filter  outputs  for  return  from  a 
range-azimuth  cell  not  containing  any  target. 


5.  SUGGESTIONS  FOR  FUTURE  STUDY 

This  paper  outlines  an  initial  study  of  a  pattern-recognition 
approach  that  has  rather  broad  applicability.  This  section  suggests 
possible  extensions  of  the  work  that  has  been  discussed  here. 

Target  identification  by  track  classification  is  essentially  a 
sequential  decision  problem.  As  experience  is  gained  with  both  sig¬ 
nature  classification  and  tracking,  the  current  procedures  will  be 
refined  so  that  targets  can  be  identified  in  a  minimum  time  from  first 
detection.  Ultimately,  the  current  signature-classification  structure 
may  evolve  into  a  sequential  classifier  where  the  various  features  are 
obtained  using  recursive  estimators.  This  approach  might  reduce  both 
the  computational  and  storage  requirements  of  a  real-time  system. 
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A  fundamental  limitation  of  track  classification  involves  a  form  of 
information  overload  that  cannot  be  circumvented  by  increased  computer 
speed  or  memory  site.  A  simple  example  will  suggest  the  nature  of  the 
problem.  Let  one  large  bird  (large  enough  to  be  detected)  be  released 
and  constrained  to  fly  within  a  grid  of  10  adjacent  range  cells  by  io 
adjacent  azimuthal  cells  (100  range-azimuth  cells  altogether).  Even 
without  problems  owing  to  range  and  azimuthal  target  spillover,  every 
cell  in  this  grid  i3  likely  to  alarm  within  this  relatively  short 
period. 

Assigning  meaningful  tracks  in  this  "high  target-density  environ¬ 
ment"  (which  could  certainly  be  worse)  is  clearly  fraught  with 
problems,  quite  independent  of  computer  capabilities.  Thus,  the  objec¬ 
tive  of  track  assessment  must  be  to  eliminate  nuisance  targets  from 
consideration  as  rapidly  as  possible  so  that  a  potentially  hostile 
target  passing  through  this  grid  of  cells  will  be  correctly  classified 
and  tracked. 

Certainly,  the  classification  algorithms  should  operate  efficiently 
with  regard  to  memory  and  computational  power  requirements  (e.g.,  cycle 
time,  parallelism).  At  least  in  theory,  the  estimated  spectral  features 
for  each  cell  within  a  track  preserve  the  essence  of  the  entire  track 
history.  Thus,  the  set  of  spectral  features  (or  possibly  only  the 
recursive  feature  statistics)  and  the  history  of  the  physical  track 
movement  (e.g.,  closing  on  a  high-value  asset)  provide  a  minimal  set  of 
data  that  completely  summarizes  the  track  history  and  yet  is  efficient 
in  both  required  memory  storage  and  computational  power. 

The  feature-estimation  algorithm  discussed  in  this  report  involves 
the  use  of  the  DFT  in  developing  the  power  spectral  density  function  of 
the  target  b&ckscatter.  An  alternative  approach  that,  in  general, 
achieves  higher  frequency  resolution  than  the  DFT  method  and  does  not 
produce  the  data-weighting  function  bias  inherent  in  this  method  is  the 
autoregressive  or  maximum  entropy  method.  In  using  this  approach,  the 
assumption  is  that  the  signal  is  generated  by  a  finite  dimensional 
process  so  that  the  signal  spectrum  can  be  expressed  by  a  rational 
polynomial  model, 12 

Because  the  target  return  appears  to  be  adequately  modelled  by  a 
finite-dimensional  system  as  discussed  in  section  2,  any  return  from  the 
signal  classes  of  interest  may  lend  itself  to  high-resolution  spectral 
analysis.  One  of  the  important  aspects  of  this  approach  is  in  obtaining 
relatively  high  resolution  from  short-length  records.  Bird  return  tends 
to  become  more  nonstationary  (in  the  statistical  sense)  as  the  record 
length  increases.  Thus,  more  reliable  classification  of  bird  targets 
may  be  possible  using  autoregressive  spectral  estimation  of  short-length 
returns. 


Proceedings  of  the  RADC  Spectrum  Estimation  Workshop,  Rome  Air 
Development  Center  (May  1978). 
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Other  areas  that  need  to  be  addressed  in  detail  during  future 
studies  include 

(1)  splitting  and  merging  tracks  and  one  track  overtaking  another 
track, 

(2)  correct  association  of  targets  for  crossing  trajectories, 

(3)  "intelligent**  centroiding  of  multiple,  adjacent-cell  detections 
involving  spatial  and  frequency-domain  assessment  of  return, 

(4)  long  dwell  time  required  by  resting  targets, 

(5)  evasive  target  tactics, 

(6)  radar  masking, 

(7)  assignment  of  track  priority, 

(8)  choice  of  an  operating  mode  in  the  event  of  data  or  computation 
overload,  and 

(9)  determination  of  optimal  tracker  memory  so  that  initial  target 

classifications  at  the  limit  of  detection  range  do  not 

adversely  bias  future  assessments. 

The  following  was  excerpted  from  an  article  by  Laveen  Kanal3  to 
provide  the  proper  perspective  on  the  admittedly  high-risk  objective  of 
reliable  target  classification  based  on  radar  backscatteri 

"For  the  most  part,  however,  technique  development  has  occurred 
without  much  feedback  from  experiments,  since  meaningful  experimentation 
in  pattern  recognition  often  requires  that  significant  resources  be 
spent  on  collection,  verification  and  handling  of  large  data  bases  . .. 

In  many  application  areas,  effective  use  of  the  data  requires  close 

interaction  with  persons  knowledgeable  about  the  processes  that  generate 
the  data.  Also  required  is  a  sustained  effort  devoted  to  the  particular 
application." 

Although  a  sustained  effort  is  still  required  to  expand  the  present 
data  base,  refine  the  current  target  classifiers,  and  explore  the  limi¬ 
tations  of  the  approach,  very  encouraging  results  have  been  obtained  to 
date.  This  program  has  had  the  benefit  of  a  relatively  large,  high- 
quality  data  base,  excellent  computer  and  graphics  facilities,  and  close 
interaction  between  experimentation  and  algorithm  development. 

Kanal ,  Patterns  in  Pattern  Recognition:  1968-1974,  IEEE  Trans . 
Inf.  Theory,  IT-20 ,  6  (November  1974),  697-722. 
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6.  SUMMARY  AND  CONCLUSIONS 


Target-classi flea ton  algorithms  have  been  developed  using  target 
modeling  in  conjunction  with  a  large  data  base  of  labeled  radar 

return.  At  the  present  time,  target  classifiers  for  walking  humans, 
vehicles,  dogs,  birds,  and  nontranslating  targets  have  been  developed 
and  tested;  results  have  been  reported  here. 

The  structural  pattern-recognition  approach  to  classifier  design 
seems  to  be  tailored  to  the  nature  of  radar  backscatter  from  the  classes 
of  moving  targets  considered.  No  classifier  training  in  the  usual  sense 
is  required,  and  thus  the  entire  data  base  can  be  used  for  testing 

classifier  performance.  Large  computing  power  is  not  required  for 
feature  selection  and  ranking.  The  method  results  in  a  very  simple 
classifier  implementation  involving  the  fast  Fourier  transform,  simple 
logic,  and  arithmetic  computations.  Additionally,  the  classifiers  can 
be  easily  refined  as  future  data  become  available. 

A  digital  simulation  of  the  radar  processor  that  develops  range- 
azimuth  cell  alarms  was  implemented  on  a  general-purpose  digital 
computer.  Based  on  this  simulation,  a  software  interface  was  developed 
between  the  detection  and  target-classification  processing.  The  current 
results  of  this  study  show  that  relatively  simple  algorithms  appear  to 
be  adequate  for  automating  target  classification  on  a  single-cell 

basis.  The  target  tracker  would  use  the  single-cell  classification 

assessment  for  each  cell  within  a  target  track  and  the  measured  feature 
statistics  in  assigning  an  overall  track  classification. 

On  a  single  range-azimuth  cell  basis,  it  appears  possible  to  provide 
reasonably  reliable  identification  of  threatening  (human  and  vehicle) 
targets  and  nonthreatening  (animal,  bird,  nontranslating)  targets.  By 
using  track  as  opposed  to  single-cell  assessment,  rather  high- 
confidence-level  classification  appears  to  be  possible  for  tracks  as 
short  as  three  or  four  range-azimuth  cells. 
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APPENDIX  A.— RADAR  EQUATION  FOR  PROPAGATION  OVER  A  PERFECTLY 
REFLECTING  GROUND  PLANE 

In  this  appendix  the  radar  equation  for  propagation  over  a  perfectly 
reflecting  ground  plane  ia  derived.  Thia  analyaia  follow  the  work  of 
Skolnik1  but  ia  adapted  to  the  oaae  of  ground  aurveillanee. 

The  preaence  of  the  ground  cauaea  a  ground -ref looted  wave  that  inter- 
cepta  the  target  with  a  different  phaae  than  the  direct  path  wave  aa 
ahown  in  figure  A-i .  The  radar  ia  aaaumed  to  be  horisontally  polariaedi 
the  ground  ia  aaaumed  to  be  a  flat,  perfectly  reflecting  aurface  with  a 
reflection  coefficient  of  -1 • 

The  difference  in  path  length  between  the  direct  wave  and  the 
reflected  wave,  6,  may  be  eaaily  determined  by  noting  two  equal! tiea  in 
the  geometrical  conatruction  of  figure  A-1 ,  i.e,. 

Direct  path  length  -  AT  «  CGT  ,  R  »  ha  ,  (A-1) 

Reflected  path  length  ■  AGT  »  BGT  . 


A 


Figure  A-1 .  Geometrical  construction  for  determining  phase 
differential  between  direct  and  ground-reflected  waves  for 
a  target  above  a  perfectly  reflecting  ground  plane. 


Vf,  Skolnik,  Introduction  to  Radar  Systems,  McGraw-Hill  (1980),  Ch 


12, 
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Thus, 


<S  =  AGY  -  AT  -  BGT  -  CGT  =  2ha  sin  £  . 


(A- 2) 


Now,  since 


GT  _  AG 


(A-3) 


R  =  AT  , 


(A-4) 


sin  £  - 


( A-5) 


it  then  follows  that 


h^  +  h. 
sin  £  »  -R  -■ 


( A-6) 


h.  +  h 
t  a 


for  small  £.  Therefore, 


6  =  2h 


ht  +  ha 


( A-7) 


Hie  phase  difference,  i}»d,  between  the  ground-reflected  and  direct  waves 
is  then  i|/d  =  2tt6/A;  the  total  phase  difference,  i|>,  between  the  two  paths 
is 


\)j  =  i);d  +  tt  . 

At  the  target,  the  resultant  power  of  two  signals  with  unit  amplitude 
but  with  phase  difference  \|>  is 


2  (l  +  cos  \ji)  =  4  si 


°2  (4)  • 


(a  a) 
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Defining  $  as  the  field  strength  at  the  target  in  the  presence  of  a 
ground  plane  relative  to  the  field  strength  at  the  target  in  free  space, 
the  power  reradiated  back  to  the  radar  via  the  direct  and  reflected 
paths  is  proportional  to 


»  16  sin4 


p-.fr  +  M]  . 


(A-9) 


The  radar  equation  is  then  scaled  by  (J4,  the  propagation  factor, 
thus: 


PR  G^2°t 
PT  ( 4-rr )  3R4 


04 


6G2X2at 
4tt)  3R4 


( A- 1 0 ) 


When  is  small,  the  radar  equation  becomes 


PR  4TTG2at 
PT  a2r® 


( A- 11) 


where 

«  -  h4(ht  +  hj4  . 

The  presence  of  the  ground  plane  causes  the  elevation  coverage  to 
break  up  into  the  lobed  pattern  given  by  equation  (A-10). 

The  first  maximum  of  the  antenna  pattern  occurs  when 


.  ha(ht  +  ha)  - 

4  - re -  -  1 


This  first  lobe  will  occur  closer  to  the  ground  if  the  antenna  height, 
h_,  is  made  larger. 
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